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ABSTRACT 

We present far-IR/sub-mm imaging and spectroscopy of 49 Ceti, an unusual circumstellar disk around a nearby 
young A1V star. The system is famous for showing the dust properties of a debris disk, but the gas properties of 
a low-mass protoplanetary disk. The data were acquired with the Herschel Space Observatory PACS and SPIRE 
instruments, largely as part of the “Gas in Protoplanetary Systems” (GASPS) Open Time Key Programme. Disk 
dust emission is detected in images at 70, 160, 250, 350, and 500 /im; 49 Cet is significantly extended in the 70 pm 
image, spatially resolving the outer dust disk for the first time. Spectra covering small wavelength ranges centered 
on eight atomic and molecular emission lines were obtained, including [O i] 63 /im and [C ii] 158 /im. The C ii line 
was detected at the 5a level — the first detection of atomic emission from the disk. No other emission lines were 
seen, despite the fact that the O i line is the brightest one observed in Herschel protoplanetary disk spectra. We 
present an estimate of the amount of circumstellar atomic gas implied by the C ii emission. The new far-IR/sub-mm 
data fills in a large gap in the previous spectral energy distribution (SED) of 49 Cet. A simple model of the new 
SED confirms the two-component structure of the disk: warm inner dust and cold outer dust that produces most of 
the observed excess. Finally, we discuss preliminary thermochemical modeling of the 49 Cet gas/dust disk and our 
attempts to match several observational results simultaneously. Although we are not yet successful in doing so, our 
investigations shed light on the evolutionary status of the 49 Cet gas, which might not be primordial gas but rather 
secondary gas coming from comets. 

Key words: circumstellar matter - Kuiper belt: general - protoplanetary disks - stars: individual (49 Ceti) 
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1. INTRODUCTION 

Protoplanetary disks begin as massive, gas-rich remnants of 

their parent interstellar molecular clouds, but rapidly evolve into 

tenuous, dusty disks within about 10 Myr. These disks, called 

debris disks, differ from protoplanetary disks in that they are 
composed of secondary material recently generated by colli- 

sions between and evaporation of asteroids and comets. Plan- 

etesimal collisions presumably give rise to rocky planets, and 

therefore the younger debris disks (< 100 Myr) provide valuable 
insight into the later stages of terrestrial planet formation. 

Debris disks are gas-poor, as evidenced by a general lack 
of sub-mm carbon monoxide emission (e.g., Zuckerman et al. 

1995). However, small amounts of gas have been found in 

several debris disks, the most famous of which is Pictoris 
(e.g., Lagrange et al. 1998; Roberge et al. 2006). This primarily 
atomic gas must also be recently produced secondary material, 
as many of the species seen have short lifetimes in optically-thin 
environments (for example, the short photoionization lifetime 
of Ci). As for the dust, the ultimate source of the gas is the 
destruction of planetesimals (see Roberge & Kamp 2011 for a 
more detailed discussion). 

The well-known 49 Ceti system consists of a bright debris 
disk surrounding a single A1V star at a Hipparcos distance 


of 59 ± 1 pc (van Leeuwen 2007). The stellar age estimated 
by comparison of the star’s position on the H-R diagram to 
theoretical stellar evolutionary tracks is either 8.9*24 Myr 
from pre-main sequence tracks or 6 1 Z 1 ^ 9 Myr from post-main 
sequence tracks (Montesinos et al. 2009). Recently, Zuckerman 
& Song (2012) identified 49 Cet as a co-moving member of the 
~40 Myr old Argus Association (Torres et al. 2008). Spectral 
energy distribution (SED) fitting suggested that the 49 Cet disk 
has two distinct components, a cold outer disk and a warmer 
inner one (e.g., Wahhaj et al. 2007). The inner dust has been 
spatially resolved in mid-IR thermal emission but no emission 
features associated with /xm-sized grains are observed in the 
Spitzer IRS mid-IR spectrum of the system (Wahhaj et al. 2007; 
Chen et al. 2006). 

49 Cet is one of two debris disks that show sub-mm CO 
emission (Zuckerman et al. 1995; Dent et al. 2005; Hughes et al. 
2008); the other is HD 21997, whose CO emission was recently 
discovered (Moor et al. 2011). The CO emission from 49 Cet led 
to questions about its evolutionary status, possibly marking it as 
a rare late-stage protoplanetary disk that has dissipated most but 
not all of its primordial gas (e.g., Hughes et al. 2008). Moderate 
gas drag at such a time may be important for establishing the final 
orbital characteristics of any terrestrial planets (e.g., Nagasawa 
et al. 2005). 
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However, the relatively low fractional infrared luminosity of 
49 Cet (Lir/L* ~ 10 -3 ; Jura et al. 1998) shows that the dust 
is optically thin, which is characteristic of a debris disk. On 
the strength of the system’s relatively advanced ~40 Myr age, 
Zuckerman & Song (2012) proposed that the observed CO gas 
is not primordial but is rather coming from a massive belt of 
colliding comet-like planetesimals. At this time, however, the 
evolutionary status of 49 Cet’s gas component is uncertain. 

In this paper, we present far-IR/sub-mm imaging and 
spectroscopy of 49 Cet obtained with the Herschel Space 
Observatory (Pilbratt et al. 2010). In Section 2, the observations 
and data reduction are described. Analysis and results appear 
in Section 3: analysis of the spatially resolved 70 /xm image in 
Section 3.1, all photometry in Section 3.2, spectroscopic analy- 
sis in Section 3.3, and simple gas mass estimates in Section 3.4. 
In Section 4, we show the new 49 Cet SED, present an im- 
proved model for the stellar spectrum, and describe our simple 
dust model fitting. Preliminary modeling of the gas component 
is discussed in Section 5. An examination of the evolutionary 
status of the 49 Cet disk material appears in Section 6. Finally, 
Section 7 summarizes our primary findings and some planned 
future work. 

2. OBSERVATIONS AND DATA REDUCTION 

49 Ceti was observed with the far-IR PACS instrument 
(Poglitsch et al. 2010), as part of the “Gas in Protoplanetary 
Systems” (GASPS) Open Time Key Programme (e.g., Mathews 
et al. 2010; Dent et al. 2013). Single PACS scan maps at 
70 /xm and 160 // in (Ohs ID 1342188485) were obtained on 2009 
December 23 using the medium scan speed (20" s -1 ). The maps 
consisted of 8 scan legs with 3' lengths and 5" cross-scan steps. 
The on-source exposure time for each map was 72 s. 

We also obtained three PACS spectroscopy observations 
in chop-nod mode. The first was a RangeSpec observation 
that included six small wavelength ranges centered on the 
following transitions: CO 73 / 2 m, H 2 O 79 /urn, CO 90 /tm, 
[Oi] 145 /ttm, [C 11 ] 158 Atm, and H 2 0 180 /xm (taken 2009 
December 22, ObsID 1342188423, 251 s on-source time per 
range). The second was a LineSpec observation that covered 
two small wavelength ranges centered on the [0 1 ] 63 Atm 
and DCO + 190 Atm transitions (taken 2009 December 22, 
ObsID 1342188424, 370 s on-source time per range). The 
third was a deeper RangeSpec observation focused on the [C 11 ] 
158 Atm line (taken 2011 July 8, ObsID 1342223790, 670 s 
on-source time). 

All PACS data were calibrated with HIPE v8.2, using pipeline 
calibration scripts (Ott 2010). Final scan maps for each wave- 
length were generated with two different pixel scales: one with 
1" pixels for use in resolved imaging analysis (Section 3.1) 
and one with the native pixel scale of the PACS detectors (372 
at 70 Atm and 674 at 160 Atm) for use in photometry analysis 
(Section 3.2). Since 49 Cet is bright at far-IR wavelengths, stan- 
dard calibration of scan maps will lead to oversubtraction of 
the background during high-pass filtering. Therefore, a region 
around the source was masked before filtering. The absolute 
calibration uncertainties for the 70 and 160 Atm scan maps are 
2.64% and 4.15%, respectively. 1 ' 

PACS is an integral field spectrometer (IFS) with a 5 x 5 
array of spectral pixels (spaxels), each 974 x 974 in size. 
All spectroscopy was obtained in ChopNod mode to remove 


13 PICC-ME-TN-037: http://herschel.esac.esa.int/twiki/pub/Public/ 
PacsCalibrationWeb/pacs_bolo_fluxcal_report_vl.pdf 


Table 1 

49 Cet Disk Geometry from the 70 /tm Image 



49 Cet 

a Boo 

49 Cet Deconvol. 

Major axis P.A. (E of N) 

-75° ± 1° 


-76.1° ±0.5° 

Major axis (FWHM) 

874 

573 

677 

Semi-major axis (HWHM) 

250 AU 


200 AU 

Minor axis (FWHM) 

671 

479 

478 

Semi-minor axis (HWHM) 

^180 AU 


$7140 AU 

Inclination 

\v 

4^ 

44* 

0 


\v 

4^ 

4^ 

0 


Note. Estimated uncertainty of an FWHM measurement is 073 =18 AU. 


telescope and sky background emission. For each spectroscopic 
observation, 25 one-dimensional spectra were extracted, one 
from each spaxel. The pipeline calibration script applied a 
flux correction to each spectrum, accounting for the expected 
fraction of the point-spread function (PSF) falling outside the 
spaxel. Spectra with two different pixel scales were produced: 
one corresponding to the native resolution of the instrument 
(2 pixels per resolution element, i.e., Nyquist sampling) and one 
with smaller pixels (oversampled with 3 pixels per resolution 
element). All analysis was performed on the Nyquist-sampled 
data, which gave results with the highest signal-to-noise ratio 
(S/N), as expected. The oversampled spectra were used only to 
check results and for display purposes. 

In addition to the PACS data, we made use of sub-mm 49 Cet 
observations taken with the SPIRE instrument (Griffin et al. 
2010). The SPIRE small scan maps at 250 Atm, 350 //m, and 
500 // m (ObsID 1342236226) were acquired in 2012 January as 
part of a Herschel Open Time program (OTl_pabraham_2, PI: 
P. Abraham). The data were calibrated using HIPE v8.2 and the 
standard pipeline calibration script, producing maps with units 
of Jy beam -1 , which we converted to Jy pixel -1 before analysis. 
The map pixel scales were 6", 10", and 14" at 250 /xm, 350 //in, 
and 500 /xm, respectively. The estimated calibration uncertainty 
for the SPIRE images is 7%. 14 

Finally, we re-calibrated the low-resolution Spitzer IRS spec- 
trum of 49 Cet published in Chen et al. (2006) using the FEPS 
Spitzer Legacy Project pipeline. As is typical, there was a flux 
disjoint between the short-low and long-low portions of the 
spectrum. We shifted the long-low portion to match the Spitzer 
MIPS 24 Aim flux (discussed in Section 3.2), which removed the 
disjoint between the two halves of the spectrum. 

3. ANALYSIS AND RESULTS 
3.1. Imaging 

49 Cet is significantly extended in the 70 /xm image compared 
to a similar image of the diskless PSF reference star a Boo 
obtained using the same medium scan speed, as shown in 
the left-hand panel of Figure 1. This is the first time thermal 
emission from dust in the outer disk of 49 Cet has been 
spatially resolved. We began our analysis by fitting two- 
dimensional (2D) Gaussians to the 49 Cet and a Boo images; 
the results appear in Table 1. We estimated the uncertainty in 
our FWHM measurements by performing the same analysis 
on a different 70 /xm scan map of a Boo. The FWHM of 
the best-fitting 2D Gaussian varied by 073 between the two a 
Boo images. Adopting this value as our FWHM uncertainty, 
49 Cet is resolved along the major axis at the 10a level. 


14 SPIRE Observer’s Manual v2.4 (2011 June): 
http://herschel.esac.esa.int/Docs/SPIRE/html/spire_om.html. 
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Figure 1. Herschel PACS 70 fim images of 49 Cet. In both panels, the pixel scale is l". Brightness contours at 3, 6, 9, and 14 times the background noise are overlaid 
(solid black contours). The dashed white ellipse shows the half-maximum contour for the best-fitting 2D Gaussian. The center of the best-fitting 2D Gaussian is marked 
with a black plus. Left: original calibrated scan map, with a 70 /im scan map of a PSF reference star ( a Boo) in the lower right corner. The position angle (PA.) of the 
49 Cet disk major axis (—75° E of N) is indicated with the straight dashed line. Right: deconvolved scan map. The new PA. of the disk major axis (—76° E of N) is 
indicated with the straight dashed line. 

(A color version of this figure is available in the online journal.) 


The measured position angle (P.A.) of the disk major axis 
(—75° ± 1° E of N) is in reasonable agreement with the P.A. 
determined from spatially resolved sub-mm CO contours (—79° 
E of N; Hughes et al. 2008). 

To estimate the true extent of the dust disk, we deconvolved 
the 49 Cet image using the a Boo image as the reference PSF, 
rotated to match the roll angle of the telescope at the time of the 
49 Cet observation. The stellar contribution was removed from 
the 49 Cet image before deconvolution, by subtracting the PSF 
reference image scaled so that the total flux from a Boo matched 
the expected stellar flux from 49 Cet. Since the expected 49 Cet 
stellar flux is only 0.2% of the actual measured flux (see 
Section 4.1), this had a negligible effect on the deconvolution. 
We employed the Richardson-Lucy deconvolution method, 
conservatively limiting the number of iterations to five. The 
reason for choosing this number will become apparent later. The 
resulting deconvolved image appears in the right-hand panel of 
Figure 1. 

The deconvolved disk image shows no sign of a central 
clearing or any asymmetric structure. The P.A. of the major 
axis changed only slightly. We re-fit a 2D Gaussian to the 
deconvolved 49 Cet image; the results are included in Table 1. 
The half-width at half-maximum (HWHM) along the major axis 
is ~200 AU. In their modeling of the 49 Cet gas and dust disk, 
Hughes et al. (2008) found that they had to set the outer radius of 
the disk to 200 AU. While the agreement between that value and 
our HWHM is partly coincidental, it does indicate that the cold 
CO gas observed in the sub-mm and the cold dust we observed 
here may be co-spatial. 

During the deconvolution procedure, we found that we could 
make 49 Cet significantly smaller than a Boo along the minor 
axis by iterating many times. But no matter how many iterations 
were used (up to ~ 1 000), 49 Cet was always larger than a Boo 
along the major axis. We took this as an indication that the 
disk is significantly resolved along the major axis but not the 
minor axis. Therefore, we chose the number of deconvolution 
iterations such that the FWHM of the disk along the minor axis 
nearly matched that of a Boo. This was accomplished in five 
iterations. Since the minor axis is not resolved, we only obtain 
a lower limit on the outer disk inclination of i ^ 44° , assuming 


the disk is inherently circular. This is broadly consistent with 
the high inclination derived from the sub-mm CO contours 
(i = 90° ± 5°; Hughes et al. 2008). 49 Cet is not significantly 
extended in the 160 /xm image or the SPIRE images. 

3.2. Photometry 

To measure the total flux from 49 Cet at 70 and 160 /xm, 
aperture photometry was performed on the scan maps with the 
native pixel scale. The radii of the circular object apertures 
were 14" and 22" at 70 and 160 /xm, respectively. These are 
larger than recommended by the Herschel PACS Instrument 
Control Centre (ICC) for aperture photometry of unresolved 
faint sources, since 49 Cet is bright and the disk is spatially 
resolved at 70 /xm. The apertures were chosen to be roughly 
1 .5 times the FWHM of the disk at that wavelength (Section 3. 1), 
so that they encompass nearly all of the disk flux. Background 
subtraction was performed using the mean brightness in sky 
annuli 20"-26" from the star center for the 70 /xm image and 
28"-34" for the 160 /xm image. Aperture corrections provided 
by the Herschel PACS ICC 1 ’ were applied to the total fluxes, 
but color corrections were not applied. 

Assuming background-limited imaging, the statistical uncer- 
tainties in the final PACS fluxes are given by 


°stat 



(1) 


where er rms is the standard deviation of the pixels in the sky 
annulus, a cor is the aperture correction, x cor is the correlated 
noise correction 1 ’ (0.95 at 70 /xm, 0.88 at 160 /xm), n. dp is the 
number of pixels in the object aperture, and « sky is the number of 
pixels in the sky annulus. The absolute calibration uncertainties 
(2.64% at 70 /xm, 4.15% at 160 /xm) were added in quadrature 
to the statistical uncertainties to give the total flux uncertainties. 
The total 49 Cet fluxes at 70 and 160 /xm are Fj o Mm = 
2. 142 Jy ± 0.058 Jy and F 160Mm = 1 .004 Jy ± 0.053 Jy. 

We performed aperture photometry on the SPIRE maps using 
circular apertures with radii of 22", 30", and 42", respectively. 
For each wavelength, the background was estimated in a sky 
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Figure 2. 49 Cet spectra showing the regions around the [Oi] 63 fim line (left) and the [Cn] 158 fim line (right). The rest wavelengths of the lines are marked with 
dotted lines. No significant O i emission is seen, while the C n line is detected at the ~5cr level. Our best-fitting models for the continuum or the continuum plus 
emission line are overplotted with red solid lines. The data shown here have the oversampled wavelength scale, highlighting the shape of the [Cn] 158 fi m emission 
line, which is unresolved. 

(A color version of this figure is available in the online journal.) 


annulus 60"-90" from the star center. Aperture corrections were 
applied and the fluxes color-corrected for a F v oc v 2 point-source 
spectrum. 15 The total 49 Cet fluxes at 250, 350, and 500 /xm are 
^250^™ = 0.372 Jy ± 0.027 Jy, Fsso^m = 0.180 Jy ± 0.014 Jy, 
and Fsoo/j.m = 0.086 Jy ± 0.009 Jy. The flux uncertainties 
include statistical and calibration uncertainties; since 49 Cet 
is relatively bright, the latter dominate. 

For use in the 49 Cet SED and to correct the Spitzer IRS 
long-low spectrum, we also calculated the continuum flux at 
24 /xm from archival Spitzer MIPS imaging taken as part of the 
MIPS-GTO program DISKLEGACY^ (PI: G. Rieke, 

AOR 21942016). The post-basic-calibrated-data image was 
downloaded from the Spitzer Heritage Archive and aperture 
photometry performed using the aperture size and correction 
from Su et al. (2006). We found F 24lim = 0.259 Jy ± 0.010 Jy, 
where the error is given by the 4% absolute calibration 
uncertainty 1 ' 1 (the statistical uncertainty is negligible). 

3.3. Spectroscopy 

In all the 49 Cet PACS spectra, significant continuum or line 
emission appears only in the central spaxel. We verified that 
the star was well-centered on the array during the spectroscopic 
observations (shifts ^07 5 ^ 0.05 spaxel). This was done by 
taking theoretical PSFs, offsetting them on a virtual PACS IFS, 
and comparing the fraction of flux in each spaxel as a function of 
shift to the observed values. All further analysis discussed here 
was performed on the spectra from the central spaxel, using the 
Nyquist-sampled data. Only one emission line is detected: [C 11 ] 
158 /x m (discussed below). 

The spectra without significant emission lines were analyzed 
in the following manner. For each spectrum, we least-squares 
fit a first-degree polynomial to the continuum. The statistical 
flux uncertainties were then estimated by taking the standard 
deviation of the fluxes minus the continuum fit, in a range 
centered on the expected line position. The width of the range 
was chosen so that roughly 68% of the pixels were within lcr of 


15 SPIRE Data Reduction Guide v2.0 (2012 March): 

hup://hcrsche!.csac.csa.int/hcss-doc-8.0/load/spirc drg/html/spirc drg.html. 

16 MIPS Instrument Handbook v3.0 (201 1 March): 
http://irsa.ipac.caltech.edu/data/SPITZER/docs/mips/ 
mipsinstrumenthandbook. 


Table 2 

Spectroscopy Results for 49 Ceti 


Line 

V 

(/xm) 

17 b 
*cont 

(Jy) 

Fi n. C 
(x 10~ 18 

a Fin , 

W m -2 ) 

S/N d 

Oi 

63.185 

2.09 ± 0.35 

<11.05 

3.68 


CO J=36-35 

72.843 

1.95 ±0.32 

<14.88 

4.96 


0 -H 2 O 

78.741 

1.90 ±0.31 

<11.86 

3.95 


CO J=29-28 

90.163 

1.88 ±0.32 

<9.25 

3.08 


Oi 

145.535 

1.16 ± 0.18 

<6.16 

2.06 


C 11 shallow 

157.741 

1.13 ±0.20 

7.23 

2.22 

3.3 or 

deep 


0.98 ±0.13 

3.690 

0.798 

4.6cr 

o-h 2 o 

179.527 

<0.73 

<12.17 

4.06 


DCO + 

189.570 

<1.37 

<13.25 

4.42 



Notes. 

a Line center rest wavelength. 

b Continuum flux at line center rest wavelength; upper limits are 3 a . The flux 
errors include statistical and absolute flux calibration uncertainties. 
c Integrated emission line flux; upper limits are 3 o'. 
d Significance of emission line detection. 

the continuum. The continuum flux values at the expected line 
center wavelengths appear in Table 2. The flux errors given are 
the statistical flux uncertainties and the absolute flux calibration 
uncertainties 1 7 ( 1 1 % at 63 /xm, 1 2% at other wavelengths) added 
in quadrature. 

Upper limits on the total emission line fluxes were calculated 
by integrating the continuum-subtracted spectrum over a small 
wavelength range centered on the expected line wavelength and 
propagating the final statistical flux errors; the wavelength range 
was ±1.5x the expected width of an unresolved emission line. 
The results appear in Table 2. A plot of the region around the 
undetected [0 1 ] 63 /xm line appears in the left-hand panel of 
Figure 2; this is the brightest line observed from protoplanetary 
disks with Herschel (e.g., Meeus et al. 2012; Dent et al. 2013). 

The only line detected, [C 11 ] 158 /xm, was seen in the original 
shallow spectrum, though it was not a significant detection 
in the data calibrated with earlier versions of HIPE (e.g., v4.2). 
Therefore, we obtained a deeper follow-up spectrum to confirm 


17 PACS Observer's Manual v2.4 (201 1 December): 
http://herschel.esac.esa.int/Docs/PACS/html/pacs_om.htmI. 
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the detection, shown in the right-hand panel of Figure 2. For both 
the shallow and deep spectra, we analyzed the [Cn] 158 /zm 
line in the following way. We modeled the continuum plus line 
by least-squares fitting the sum of a first-degree polynomial 
and a Gaussian. Then the statistical flux uncertainties were 
estimated by taking the standard deviation of the fluxes minus 
the best-fitting total model, in a range around the line. This time, 
the width of the range was chosen to maximize the S /N of the 
emission line (±7.5 x the instrumental FWFIM). As before, the 
total C ii emission line flux was integrated from the continuum- 
subtracted data; the results appear in Table 2. The integrated 
flux from the deep spectrum is within 3 a of the flux from the 
noisier shallow spectrum. 

To date, 49 Cet is one of only two disks observed with 
Herschel that show Cii emission but no Oi. The other is 
HD 32297 , a bright ~30 Myr old debris disk system (Donaldson 
et al. 2013; Kalas et al. 2005). HD 32297 was previously known 
to contain atomic gas, but no molecular gas has yet been detected 
(Redfield 2007). 


3.4. Carbon Gas Mass 

With only one carbon emission line detected, we are not able 
to measure the excitation temperature of the gas from our data, 
which is needed to calculate a model-independent total mass of 
Cii ions from the integrated [Cii] 158 /z m line flux. However, 
we may determine a temperature-insensitive lower limit on the 
total mass, assuming the emission is optically thin. In this case, 
the mass is given by 


Men = 


Arc Xq 
he 


Tint m d 2 

A ui %u 


( 2 ) 


where Xq is the wavelength of the line, /± lt is the observed 
integrated emission line flux, m is the mass of an atom, cl is 
the distance from the emitting region to the observer, u and 
/ designate the upper and lower energy levels involved in the 
transition, A u i is the spontaneous transition probability, 18 and 
x u is the fraction of atoms in the upper energy level. Assuming 
local thermal equilibrium (LTE), x u is given by 

(2 J u + 1) e- £ «/* r “ 


where J u is the angular momentum quantum number of the upper 
level, E u is the energy of the upper level, 7' cx is the excitation 
temperature, and Q Ta is the partition function for the given 
excitation temperature. 

Figure 3 shows a plot of the C ii mass as a function of the 
assumed excitation temperature. The lower limit on the mass 
is Men > 2.15 x 1 0 4 M®, valid for T ex between 1 K and 
2000 K. This value is close to the total CO mass calculated from 
the observed sub-mm CO emission ( Mqq = 2.2 x 10~ 4 M®; 
Hughes et al. 2008). Our upper limits on the CO 73 and 90 fim 
line fluxes do not provide useful checks on the CO abundance or 
excitation temperature. These lines arise from very high energy 
levels ( E u = 3471 and 2240 K, respectively) and therefore are 
insensitive tracers of the bulk of the CO gas. 

To determine the total mass of carbon atoms in the 49 Cet disk, 
we need to know the ionization balance in the gas. This may be 
calculated with complex thermochemical disk models, further 
discussed in Section 5. Here we make a simple estimate of the 


18 NIST Atomic Spectra Database: http://www.nist.gov/pml/data/asd.cfm. 
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Figure 3. 49 Cet C ii gas mass inferred from the 158 /cm emission line flux as a 
function of excitation temperature. The ± ler errors on the mass are plotted with 
green lines. The horizontal dotted line shows the lower limit on the gas mass 
(Men > 2.15 x ! 0 4 M®), valid for assumed excitation temperatures between 
1 K and 2000 K. 

(A color version of this figure is available in the online journal.) 


mass lower limit by assuming a plausible ionization fraction, 
measured in an analogous environment. The radiation and 
density characteristics of the 49 Cet circumstellar environment 
are similar to those of fi Pic; both A stars are surrounded by 
optically thin disks containing a similar amount of dust, judging 
from the systems’ fractional infrared luminosities. Both disks 
contain circumstellar gas, including CO (although the /I Pic CO 
is not yet detected in sub-mm emission, only UV absorption; 
Roberge et al. 2000). 

The measured column densities of neutral and first-ionized 
carbon in the Pic gas are roughly equal, giving a carbon ioniza- 
tion fraction of ~50% (Roberge et al. 2006), in agreement with 
the value predicted by photoionization calculations (Fernandez 
et al. 2006). Therefore, the rough lower limit on the total mass 
of carbon atoms in the 49 Cet gas is Me > 4.3 x 10 -4 AT®. 
If we further assume that the atomic gas has solar elemental 
abundances, we may estimate a lower limit on the total atomic 
gas mass by dividing the carbon mass limit by the solar carbon 
mass fraction (0.288%; Lodders 2003). With all these assump- 
tions (optically thin [C ii] emission, LTE, carbon 7' cx < 2000 K, 
50% ionization fraction, and solar elemental abundances), the 
rough lower limit on the total atomic gas mass is M a >0.15 M® . 

Turning to the other undetected emission lines, we cannot at 
present use them to calculate model-independent upper limits on 
the masses of other gases, due to a lack of precise information on 
the gas excitation temperatures. Examination of Figure 3 shows 
that the C ii mass versus excitation temperature curve changes 
slowly at high temperatures but shoots up at low temperatures; 
this is a generic characteristic of such curves. It is easy to set a 
mass lower limit that is valid for a wide range of temperatures, 
since it is controlled by the highest temperature considered, 
where the curve is changing slowly. Setting a temperature- 
insensitive upper limit is difficult, since it is controlled by the 
lowest temperature considered, where the curve is changing 
rapidly. 

The situation is worse for the [O i] 63 /im line, since the 
energy of the upper level ( E u = 228 K) is much higher 
than it is for the [Cii] 158 /im line ( E u = 91 K). Therefore, 
the Oi mass upper limit is extremely sensitive to the range 
of temperatures considered. For example, assuming LTE and 
changing the lowest excitation temperature considered from 
10 K to 50 K decreases the Oi mass upper limit by six orders 
of magnitude. Furthermore, the large energy separation for the 
levels involved in the [O i] 63 jim line means that the levels are 
unlikely to be thermalized and assuming LTE is not appropriate. 
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Table 3 

Unresolved Continuum Fluxes for 49 Cet 


Source 

Band 
or Mode 

Wave. 

(pm) 

Flux 

(Jy) 

Reference 

Tycho-2 

Johnson B 

0.44 

22.26 ± 0.35 

Hog et al. (2000) 

Tycho-2 

Johnson V 

0.55 

21.60 ±0.20 

Hog et al. (2000) 

2MASS 

j 

1.24 

10.18 ±0.19 

2MASS All-Sky Catalog of Point Sources 

2MASS 

H 

1.65 

6.30 ±0.13 

2MASS All-Sky Catalog of Point Sources 

2MASS 

K, 

2.16 

4.373 ± 0.081 

2MASS All-Sky Catalog of Point Sources 

WISE a 

Wl 

3.35 

1.99 ±0.12 

Wright et al. (2010) 

WISE a 

W2 

4.60 

1.294 ± 0.049 

Wright et al. (2010) 

WISE a 

W3 

11.56 

0.211 ±0.021 

Wright et al. (2010) 

IRAS b 


12.00 

0.228 ± 0.034 

IRAS Faint Source Catalog, v2.0 

Keck/MIRLIN 


12.50 

0.200 ± 0.026 

Wahhaj et al. (2007) 

Keck/MIRLIN 


17.90 

0.186 ±0.025 

Wahhaj et al. (2007) 

WISE a 

W4 

22.09 

0.238 ± 0.024 

Wright etal. (2010) 

Spitzer/MlPS c 


24.00 

0.259 ± 0.010 

This work 

IRAS b 


25.00 

0.312 ±0.042 

IRAS Faint Source Catalog, v2.0 

IRAS b 


60.00 

2.17 ±0.12 

IRAS Faint Source Catalog, v2.0 

IRAS b 


100.00 

1.88 ±0.21 

IRAS Faint Source Catalog, v2.0 

Herschel/ PACS d 

Spec. 

63.19 

2.09 ± 0.35 

This work 

Herschel/ PACS d 

Phot. 

70.00 

2.142 ±0.058 

This work 

Herschel/ PACS d 

Spec. 

72.84 

1.95 ±0.32 

This work 

Herschel/ PACS d 

Spec. 

78.74 

1.90 ±0.31 

This work 

Herschel/ PACS d 

Spec. 

90.16 

1.88 ±0.32 

This work 

Herschel/ PACS d 

Spec. 

145.54 

1.16 ±0.18 

This work 

Herschel /PACS d 

Spec. 

157.68 

0.98 ±0.13 

This work 

Herschel /PACS d 

Phot. 

160.00 

1.004 ± 0.053 

This work 

Herschel/ PACS d 

Spec. 

179.53 

<0.73 (3<t) 

This work 

Herschel/ PACS d 

Spec. 

189.57 

<1.37 (3<t) 

This work 

Herschel / SPIRE 1 * 

Phot. 

250.00 

0.372 ± 0.027 

This work 

Herschel/ SPIRE d 

Phot. 

350.00 

0.180 ±0.014 

This work 

Herschel / SPIRE d 

Phot. 

500.00 

0.086 ± 0.009 

This work 

JCMT 


800.00 

<0.039 (3a ) 

Zuckerman & Becklin (1993) 

JCMT/SCUBA 


850.00 

0.0082 ± 0.0019 

Song et al. (2004) 

IRAM 


1200.00 

0.0127 ± 0.0028 

Bockelee-Morvan et al. (1994) 

SMA 


1300.00 

<0.01 (3(7 ) e 

Hughes et al. (2008) 


Notes. 

a WISE Wl, W2, and W3 fluxes are color-corrected for a v 2 spectrum. The W4 flux is color-corrected for a 
100 K blackbody. The WISE errors include statistical and absolute calibration uncertainties added in quadrature (see 
http://wise2.ipac.caltech.edu/docs/release/allsky/expsup/sec4__4h.html). 

b IRAS fluxes are color-corrected (see http://lambda.gsfc.nasa.gov/product/iras/colorcorr.cfm). 
c For the Spitzer MIPS 24 pm flux, the error is the 4% absolute calibration uncertainty. 

d For all Herschel data, the flux errors are the statistical and calibration uncertainties added in quadrature. The SPIRE 
fluxes are color-corrected for a v 2 spectrum (see http://herschel.esac.esa.int/Docs/SPIRE/html/spire_om.html). 
e Upper limit on total continuum flux, assuming the dust emission is spread over 4.8 SMA beams (beam size = l'.'O x 172). 
This is the size of the CO 7 = 2—1 emission region in the SMA images (Hughes et al. 2008). 


4. SED FITTING 

We have collected a wide range of unresolved photometry 
and spectroscopy of 49 Cet; our compilation of continuum 
fluxes at wavelengths from the optical to the millimeter appears 
in Table 3. In Figure 4, we show the most complete SED to 
date. It includes our new far-IR/sub-mm continuum fluxes from 
Herschel, the Spitzer IRS spectrum; the Spitzer MIPS 24 /tm 
flux; and new near-IR photometry from WISE. Some redundant 
fluxes given in Table 3 were not used in the SED. 

Table 3 includes two unrefereed detections of 49 Cet at sub- 
mm/mm wavelengths: one at 850 pm from the lames Clerk 
Maxwell telescope (JCMT)/SCUBA (Song et al. 2004) and one 
at 1.2 mm from IRAM (Bockelee-Morvan et al. 1994). As noted 
by Hughes et al. (2008), the two fluxes cannot be simultaneously 
fitted with any dust excess model that decreases with increasing 
wavelength in the sub-mm/mm. It is not obvious which flux 
value is more accurate. In the absence of other sub-mm fluxes, 


Hughes et al. (2008) conservatively adopted the lower flux value 
from SCUBA (Song et al. 2004). The new SPIRE fluxes allow 
us to reasonably constrain the SED fitting without relying on the 
ambiguous SCUBA and IRAM fluxes. Therefore, we avoided 
making an arbitrary choice between them by excluding both 
points from our new SED. We eagerly await sensitive new 
sub-mm/mm observations of 49 Cet with ALMA. 

4.1. Stellar Model 

To determine the expected stellar flux at all wavelengths, we 
compared photospheric models to the 49 Cet SED shortward of 
about 8 /xm, where the dust excess starts to become noticeable. 
A careful analysis of the 49 Cet stellar parameters appears in 
Montesinos et al. (2009). We were motivated to revisit that 
work by the existence of important information not available 
or utilized at the time: specifically, the independent stellar age 
estimate (~40 Myr; Zuckerman & Song 2012) and the Spitzer 
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Wavelength (pm) 

Figure 4. Spectral energy distribution for 49 Cet with a two-component 
blackbody model for the dust excess emission. Our best-fitting ATLAS9 stellar 
model (solid black line) was fit to the B and V magnitudes from the Tycho-2 
Catalog, the J, H, and K s fluxes from 2MASS, and the WISE W1 flux (cyan 
circles). The total model (black dashed line) is the sum of the best-fitting stellar 
model, a warm simple blackbody, and a cold modified blackbody. It was fit to 
various photometry detections (red diamonds), continuum fluxes from PACS 
spectroscopy (gray diamonds), and the Spitzer IRS spectrum (green line). The 
individual blackbody components are plotted with orange dotted lines. Some 
relevant flux upper limits are also shown (blue arrows). 

IRS spectrum. We found that while the photospheric model 
computed using the best stellar parameters in Montesinos et al. 
(2009) fit the optical and near-IR fluxes well, it appeared slightly 
too bright when compared with the short wavelength portion of 
the IRS spectrum. 

Therefore, we fit ATLAS9 solar-metallicity photospheric 
models 19 to the optical and near-IR data (Castelli & Kurucz 
2004). The data used were the B- and l'- band fluxes from the 
Tycho-2 Catalog, the J, H, and K s fluxes from the Two Micron 
All Sky Survey (2MASS), and the color-corrected WISE Band 1 
flux. No interstellar reddening was applied to the models. The 
best-fitting model had an effective temperature T e ff = 9000 K 
and also matched the short wavelength portion of the IRS 
spectrum well. This temperature is typical of an A1V star 
but slightly cooler than the 9500 K effective temperature in 
Montesinos et al. (2009). However, this fitting procedure could 
not provide a useful constraint on the stellar gravity. 

Using the new effective temperature and age estimate, we 
placed the star on a log g— log T e g H-R diagram, forcing it to lie 
precisely on the 40 Myr isochrone. We used the evolutionary 
tracks and isochrones with solar metallicity from the Yonsei- 
Yale group (Yi et al. 2001). The stellar gravity indicated by the 
age constraint is log g = 4.3, the corresponding stellar mass 
is 2 M 0 , and the stellar luminosity is 15.5 L Q . A photospheric 
model with T e ^ — 9000 K and log g = 4.3 was built by interpo- 
lating two ATLAS9 solar-metallicity photospheric models with 
log g — 4.0 and 4.5. No extinction was applied to the model, 
which was then normalized to the IRS spectrum at 6 pt m. The 
final stellar model is shown with a black solid line in Figure 4. 

A strong check of the internal consistency of the set of stellar 
parameters consists of using the stellar luminosity derived from 
the age constraint (L* = 15.5L 0 ) and the total photospheric 
flux (F*, computed by integrating the final stellar model) to 


19 Downloaded from F. Castelli’s public database at 
http://wwwuser.oat.ts.astro.it/castelli/grids.html. 


Table 4 

Parameters of the Best Blackbody Dust Models 


Parameters 

S ingle-temperature 

Two-temperature 

T (K) 

83 ± 9 

62 ±1 175 ±3 

ko (/im) 

156 ±35 

102 ±12 

P 

0.6 ±0.1 

0.7 ±0.1 

£dus./T* (X 10~ 4 ) 

11 ± 2 

7.1 ±0.2 1.26 ±0.03 

Reduced y 2 

4 x 10 6 

1.18 


recover the observed distance. Using the simple expression 
L* = 47r d 2 , the derived distance is 56 pc, in excellent 
agreement with the Hipparcos distance of 59 pc. A lower gravity 
would imply a larger luminosity and hence a much greater 
distance; as an example, Teff = 9000 K and log g = 4.0 in the 
log L*/L 0 — log T e a H-R diagram corresponds to a luminosity 
of 34.67 L q (using PMS tracks) or 36.31 L 0 (using post- 
MS tracks). These luminosities give distances of 83-85 pc, 
inconsistent with the Hipparcos distance. 

4.2. Modified Blackbody Dust Models 

We began our analysis of the 49 Cet dust excess emission 
with a very simple model, i.e., the sum of the best stellar 
model and a single-temperature modified blackbody. Starting 
from Equation (6) in Beckwith et al. (1990), an optically thin 
modified blackbody has the form: 

F V (T) = B V (T) (1 - e~ z «) £2 ~ B V (T) r v £2 , (4) 

where B V (T) is the Planck function for temperature T, r v is 
the optical depth, and £2 is the solid angle of the emitting area 
(effectively a constant of proportionality). The optical depth 
is proportional to the dust mass opacity. At mid-IR to mm 
wavelengths, mass opacities for expected debris disk grains 
have a power-law dependence on wavelength (see Figure 7 in 
Lebreton et al. 2012). At shorter wavelengths, the mass opacity 
for simple astronomical silicate grains abruptly becomes nearly 
independent of wavelength. Therefore, our formulation for the 
optical depth is 

T v = To A. < 2-0 

= T 0 (i) x > x o, (5) 

where to is the optical depth at reference wavelength Xq and p 
is the emissivity power-law index. 

This simple model was least-squares fit to the 49 Cet SED, 
ignoring flux upper limits. We found that a single-temperature 
blackbody model could not simultaneously fit the mid-IR 
emission seen with Spitzer and the far-IR emission peak. The 
parameters for the best dust model that fits the far-IR data 
appear in Table 4. The fractional infrared luminosity (idust/E*), 
which is the ratio of the blackbody model integrated over 
all frequencies to the similarly integrated stellar model, is 
proportional to Tq. The reduced / 2 value is extremely large 
due to the model’s failure to also match the numerous points in 
the IRS spectrum. That failure indicates that the disk has dust 
grains at a wide range of temperatures, which can be interpreted 
as a warm inner dust component and cold outer material. Wahhaj 
et al. (2007) came to the same conclusion and fit the 49 Cet SED 
with a two-component dust disk. 

We therefore made a total SED model consisting of the best 
stellar model, a warm simple blackbody, and a cold modified 
blackbody. This model was also least-squares fit to the 49 Cet 
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SED, ignoring flux upper limits as before. The best-fitting two- 
temperature model is overplotted in Figure 4, along with the two 
individual blackbody components; the model parameters appear 
in Table 4. The reduced y 1 value is greatly improved compared 
to the single-temperature model. 

The temperature of the warm component (175 ± 3 K) is 
higher than the temperature found by fitting a single-temperature 
simple blackbody to the IRS spectrum alone (118 ± 6 K; 
Chen et al. 2006). This is presumably because we are now 
accounting for the flux in the IRS spectrum coming from the 
cold component. The emissivity power-law index found for the 
cold component = 0.7 ± 0.1) is smaller that the values 
typically seen in the interstellar medium ~ 2; e.g., Andrews 
& Williams 2005). This is often seen in protoplanetary and 
debris disk SEDs and it indicates that the emitting grains are 
larger than the typical sub-micron interstellar grains (e.g., Su 
et al. 2006). 

Assuming the dust grains are in radiative equilibrium with the 
central star, we estimated the minimum characteristic radius of 
each dust component. Starting from Equation (1) in Beckwith 
et al. (1990), we adopted simple blackbody grains and used the 
bolometric stellar flux from our best stellar model. More realistic 
grains will be further from the star than simple blackbody grains 
at the same temperature, so the distances calculated are lower 
limits on the characteristic dust radii. The characteristic radius 
of the warm component is >11 AU, while that of the cold 
component is >84 AU. 

5. GAS MODELING 

There are two possible explanations for the gas content of the 
49 Cet disk. As mentioned in the Introduction, the most common 
previous scenario interprets the CO as the last remnant of the 
primordial gas left over from star formation (e.g., Hughes et al. 
2008). In the first subsection, we describe our attempts to fit the 
49 Cet SED and gas observations using a primordial-like disk 
model. The next subsection considers the second explanation, 
that the 49 Cet gas is secondary material coming from the 
destruction of planetary material. 

5.1. Primordial Gas Scenario 

We began a detailed investigation of the 49 Cet gas emission 
using an advanced disk modeling code, ProDiMo, which self- 
consistently determines the temperature, density, and chemical 
structure of combined gas and dust disks (Woitke et al. 2009; 
Kamp et al. 2010; Thi et al. 2011). ProDiMo calculates the 
heating and cooling of both gas and dust, and takes into account 
960 chemical reactions involving 71 species. Photochemistry, 
H 2 formation on grain surfaces, and cosmic-ray heating are all 
included. 

As a first step, we considered the 49 Cet disk model from 
Hughes et al. (2008) that successfully reproduced the basic 
features of the spatially and spectrally resolved sub-mm CO 
emission, as well as the dust emission from the outer disk (the 
inner disk was ignored). In that work, steady-state chemistry 
was adopted and solar elemental abundances assumed. Some 
key features of the model were (1) a lack of gas within 40 AU 
of the star, (2) an outer disk radius of 200 AU, and (3) a total 
gas mass of 13 M®. 

In the original Hughes et al. (2008) model, the gas-to-dust 
ratio had to be set to the unusually high value of 500. This 
was likely caused by necessary simplifying assumptions about 
the dust properties: well-mixed gas and dust, a uniform dust 


composition (astronomical silicates), and a single large grain 
size. Here we retained the first two assumptions but used a 
Dohnanyi power-law grain size distribution ( dn oc s~ 35 ds; 
Dohnanyi 1969). We then calculated several models that re- 
produced the observed sub-mm CO emission while varying 
other disk parameters (e.g., gas-to-dust ratios from 300 to 
65). While it was possible to find models that matched the 
sub-mm CO emission and the upper limit on [0 1 ] 63 /r m emis- 
sion, they underpredicted the [C 11 ] 158 /u m emission by factors 
of 5-23. 

In sum, simple adjustments to the Hughes et al. (2008) model 
were unsuccessful at simultaneously reproducing the sub-mm 
CO emission, the Cii emission, and the lack of Oi emission. 
For example, reducing the oxygen abundance in an attempt to 
decrease the 0 1 emission led to less CO formation and reduced 
CO cooling, which resulted in a warmer gas disk. Paradoxically, 
this had the effect of increasing the Oi emission rather than 
decreasing it. The only adjustment that showed some promise 
involved increasing the carbon abundance over the solar value. 
Further work on detailed modeling of the 49 Cet disk with 
ProDiMo is underway. For now, it appears that the observations 
cannot be reproduced with a “primordial-like” protoplanetary 
disk model. 


5.2. Secondary Gas Scenario 

We therefore consider whether the gas could also be con- 
sistent with a secondary source from planetesimals. If so, its 
composition could shed light on the make-up of the young par- 
ent planetesimals — outcomes of the planet formation process 
and the building blocks of full-sized planets. Questions that 
can be addressed include whether the parent planetesimals are 
rocky or icy and whether there are any abundance anomalies 
(e.g., Roberge et al. 2006; Xie et al. 2013). 

In this scenario, the presence of relatively abundant CO (for a 
debris disk) suggests that this gas comes from icy material, more 
comet-like than asteroid-like (as proposed in Zuckerman & Song 
2012). There are several possible mechanisms for its production, 
including outgassing of comet-like bodies, photodesorption of 
ice-coated grains, and grain-grain collisions (Lagrange et al. 
1998; Chen et al. 2007; Czechowski & Mann 2007). How large 
a mass of comets would be needed to supply the observed 
mass of CO? In equilibrium, the CO production rate will 
match the loss rate. Assuming the CO is in a low-density 
environment, chemical reactions can be ignored and the primary 
loss mechanism is photodissociation. 

Therefore, the CO loss rate is 


/ dnc o 
V dt 


= - k n C o = -k 

loss 



(6) 


where nco is the CO number volume density, k is the photodis- 
sociation rate, Xc® is the total number of CO molecules, and 
V is the total volume of the gas. Given the low Ld us t/Z* value, 
shielding of the CO from dissociating UV radiation by dust 
grains should be negligible. Shielding by H 2 is hard to assess, al- 
though most hydrogen in comets is locked in H 2 O rather than H 2 
(e.g., Mumma & Charnley 201 1). For now, we ignore the com- 
peting effects of shielding and the stellar UV field, and consider 
only the interstellar field. The unshielded photodissociation rate 
for CO in a Draine interstellar UV field is k = 2.6 x HU 1(1 s -1 , 
giving an unshielded CO lifetime Tqo = I / A: ~ 120 yr (Visser 
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dn 


CO 


dt 


/ prod. 


<^h 2 o r co 

+ /'ro — H 

dt V 


d Xh 2 o 
dt 


(7) 


where rco is the fractional abundance of CO relative to water 
in comets, dn Hl0 /dt is the water production rate, and X H2 q 
is the total number of FbO molecules. Observations of solar 
system comets show rco values ranging from about 0.4% to 
30% (Mumma & Charnley 201 1). 

Equating the production and loss rates shows that equilibrium 
is achieved for a water mass loss rate of 

d_M^ = J<_ iniho Mcq ' (8) 

dt r C o tn co 

where Mh 2 o is the total mass of H 2 O, /wh 2 o is the mass of a hCO 
molecule, mco is the mass of a CO molecule, and Mco is the 
total mass of CO (2.2 x 1 0 4 M e ; Hughes et al. 2008). Therefore, 
the water mass loss rate needed to produce the observed CO is 


dMli2 ° sa i(T 13 to 1CT 11 M e s“‘ . (9) 

dt 

Since comets are mostly water, we can adopt these values as 
approximate comet mass loss rates. Collisional modeling of the 
solar system’s Kuiper Belt predicts an initial mass of ~60 M e 
(Campo Bagatin & Benavidez 2012). Scaling this value up 
by the larger stellar mass of 49 Cet («»2 M 0 ) gives an initial 
planetesimal mass of ~ 120 M e . With the loss rates given above, 
it would take between roughly 0.4 Myr and 32 Myr to exhaust 
the total initial planetesimal mass. Since the star is likely to be 
~40 Myr old (Zuckerman & Song 2012), this analysis would 
suggest that if the CO is being produced from comets, it is a 
relatively short-lived phenomenon. 


6. DISCUSSION 

The dissipation of abundant primordial material left over 
from star formation sets crucial constraints on the formation 
of planetary systems. On one hand, the removal of the gas limits 
the time available for formation of gas giant planets. On the 
other hand, the presence of a modest amount of gas during the 
later stages of formation can help damp the inclinations and 
eccentricities of terrestrial planets (e.g., Nagasawa et al. 2005). 
Therefore, understanding the speed of gas dissipation and the 
mechanisms by which it occurs is vitally important for informing 
planet formation theories. 

Primordial dust lifetimes are fairly well-constrained by ob- 
servations to be <10 Myr, although there is a large dispersion 
in values for individual stars of nearly the same age and uncer- 
tainties remain about the effects of stellar mass, binarity, and 
star-forming environment (e.g., Haisch et al. 2001; Andrews 
& Williams 2005; Carpenter et al. 2006; Cieza et al. 2009; 
Luhman et al. 2010). Determining gas lifetimes, one of the pri- 
mary goals of the Herschel GASPS project, is far more difficult. 
While arguments can be made that gas and dust dissipate on 
roughly similar timescales (Roberge & Kamp 2011), there is 
significant uncertainty about the co-evolution of gas and dust 
(e.g., Pascucci et al. 2009). 

From this perspective, 49 Cet might be a nearly unique disk 
system that is just finishing dissipation of its primordial gas. 
With gas emission reminiscent of a low-mass protoplanetary 
disk, 49 Cet does look to be in some sort of intermediate state. 


As pointed out in Hughes et al. (2008), the presence of an 
inner dust disk that lacks significant molecular gas could be 
consistent with disk dissipation from the inside out through 
photoevaporation, if the inner disk grains are large enough not 
to be entrained in a photoevaporative flow (e.g., Alexander 
et al. 2006). With the previous age estimate of ~9 Myr 
(Montesinos et al. 2009), 49 Cet would have a relatively long 
but not implausible primordial gas lifetime. However, the new 
40 Myr age (Zuckerman & Song 2012) requires an anomalously 
long lifetime for the observed CO to be primordial gas. Our 
preliminary disk modeling with ProDiMo casts further doubt on 
the 49 Cet gas being primordial, since we have great difficulty 
fitting all the gas observations with a low-mass but otherwise 
normal protoplanetary disk model. 

Looking at 49 Cet’s dust properties, it greatly resembles a 
young debris disk. The low fractional dust luminosity (<1% of 
a typical Herbig Ae disk value; Meeus et al. 2012) and the lack 
of 10 jim silicate emission indicating few small grains in the 
warm inner disk are both highly characteristic of debris disks 
(e.g., Chen et al. 2006). All of the estimated ages for 49 Cet 
could place it in a late stage of terrestrial planet formation (e.g., 
Kenyon & Bromley 2006). Furthermore, recent work has shown 
that many debris disk SEDs are best fit by two-temperature 
models, like 49 Cet, suggesting radially separated inner and 
outer dust belts (Morales et al. 201 1). 

The 49 Cet gas may also be secondary debris material coming 
from destruction of planetesimals. Our simple calculations of 
gas production from comet-like ices indicate that the amount of 
material required to produce the observed CO is not implausibly 
large. Unfortunately, we are not able to say at this time if the 
lower limit on the C 11 mass and non-detection of 0 1 emission 
are also consistent with the comet-evaporation scenario. First, 
if the gas is coming from planetesimals, then the parent 
species for the bulk of the carbon gas is not necessarily CO. 
In the better- studied case of fi Pic, balancing Ci production 
by photodissociation of CO with Ci loss by photoionization 
suggests that only about 2% of the total carbon gas comes 
from dissociation of CO (Roberge et al. 2000, 2006). Since 
the primary atomic gas production mechanism is likely to be 
photodesorption of grain surfaces or grain-grain collisions, 
many other materials could supply carbon gas to the disk (e.g., 
amorphous carbon). This makes interpreting the lower limit on 
the 49 Cet C 11 mass difficult. Second, while we expect oxygen 
gas as an end-product of water ice evaporation, we are not able 
to turn the observed 0 1 flux upper limit into a useful oxygen 
mass upper limit (see Section 3.4). 

However, recent monitoring of Ca 11 lines in optical spectra of 
49 Cet show absorption features that are variable in both strength 
and velocity shift (Montgomery & Welsh 2012). Similar features 
are also seen in spectra of Pic and are attributed to so-called 
falling evaporating bodies, which are star-grazing planetesimals 
passing through the line of sight to the central star (e.g., Beust 
et al. 1990). Therefore, each absorption feature is effectively a 
transit of an exo-comet or exo-asteroid. Taken together, the 
difficulty fitting the 49 Cet observations with a primordial 
disk model, the possibility of producing the observed CO 
from comet-like material, and the detection of gas-producing 
planetesimals all indicate that 49 Cet is a close analog of 
the better-understood /I Pic debris disk. Tighter constraints on 
the gas-to-dust ratios in both systems, through more accurate 
determination of dust masses and measurements of additional 
gas species, would help confirm their similarity and allow 
studies of the effect of the gas on grain dynamics. 
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7. CONCLUDING REMARKS 

The Herschel far-IR imaging and spectroscopy of 49 Cet has 
provided an extremely rich data set for better characterizing 
both the dust and gas. Our primary findings are as follows. 

1. In our 70 fim image, dust continuum emission from the 
outer disk is spatially resolved for the first time. The disk 
orientation agrees with that determined from spatially and 
spectrally resolved sub-mm CO emission, suggesting the 
gas and dust are co-spatial. The deconvolved image does 
not show a central clearing in the dust; the Gaussian HWHM 
along the disk major axis is ~200 AU. 

2. The [C ii] 158 // m emission line was detected at the 3 a level 
in a preliminary shallow spectrum and at the 5a level in a 
deep follow-up spectrum. No other emission lines were 
detected, including [O i] 63 /xm. The integrated C ii line 
flux is (3.690 ± 0.798) x 10~ 18 W m -2 , giving an ionized 
carbon gas mass Men > 2.15 x 10~ 4 M®. With various 
assumptions, this suggests a lower limit on the total atomic 
gas mass of M a >0.15 M®. 

3. Simple modeling of the new SED confirms the two- 
component structure of the disk. We find blackbody tem- 
peratures of 175 K for the inner dust and 62 K for the outer 
dust, implying characteristic radii of >1 1 AU and >84 AU. 

4. Preliminary thermochemical modeling of the combined 
gas and dust disk indicates that a “primordial-like” pro- 
toplanetary disk model cannot simultaneously reproduce 
all the gas observations. The gas may instead be sec- 
ondary material coming from the destruction of comet-like 
ices. 

The next step in our work on the 49 Cet system is fitting 
the SED with an advanced dust disk modeling code (GRaTeR; 
Augereau et al. 1999). This code includes realistic grain proper- 
ties and can take into account the powerful constraints pro- 
vided by resolved imaging, as was done in our modeling 
of the HD 181327 debris disk (Lebreton et al. 2012). The 
new dust disk model will provide a more accurate dust mass 
and can be input to ProDiMo for further exploration of the 
disk parameter space. We will also experiment with ProDiMo 
runs in “debris disk mode,” adopting time-dependent chem- 
istry and non-solar starting compositions. Through these and 
other studies, we hope to resolve the nature of the interesting 
49 Cet disk. 

Herschel is an ESA space observatory with science instru- 
ments provided by European-led Principal Investigator consor- 
tia and with important participation from NASA. Support for 
this work was provided by the NASA Herschel Science Center 
through an award issued by JPL/Caltech. A. Roberge also ac- 
knowledges support by the Goddard Center for Astrobiology, 
part of the NASA Astrobiology Institute. J.-C. Augereau thanks 
the CNES-PNP for financial support. C. Eiroa, G. Meeus, and 
B. Montesinos were partly supported by Spanish grant AYA 
2011-26202. 
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HIPPARCOS, CTIO:2MASS, WISE, IRAS, JCMT, SMA 
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